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Abstract 

We determine the kind of parametric oscillators that are generated in the usual factorization procedure of 
second-order linear differential equations when one introduces a constant shift of the Riccati solution of 
the classical harmonic oscillator. The mathematical results show that some of these oscillators could be of 
physical nature. We give the solutions of the obtained second-order differential equations and the values of 
the shift parameter providing strictly periodic and antiperiodic solutions. We also notice that this simple 
problem presents parity-time (PT) symmetry. Possible applications are mentioned. 
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1. Riccati nonlinear equations have many applications in physics [l-3] and their solutions 
under the name of superpotentials play an important role in supersymmetric quantum me- 



chanics (4], which mathematically speaking is the result of a particular case of the factor- 
ization algorithm as applied to second-order linear differential equations. However, the 
latter equations form the groundages not only of quantum mechanics but also of many other 
branches of physics. For example, classical harmonic oscillator equations are amongst the 
simplest such equations in classical mechanics. Some time ago, factorization techniques have 
been used by Rosu and Reyes for the damped harmonic classical oscillator In addition, 



Rosu et al [7| used a Darboux-embedding technique |8[ to enclose N underdamped classical 
modes in a chirp frequency with a hyperbolic secant time profile. The main goal of the 
present letter is to investigate simple factorization procedures for the free harmonic classical 
oscillator. As it is well known, particular solutions R of the Riccati equation in its standard 
form 

R' + R 2 + f = , (1) 

where / is a function of the independent variable, enter the factorization brackets of linear 
second order differential equations 

1 + R ) ( -T - R \ u = = u " ~ ( R2 + R ') u = = u" + fu = . (2) 



dt J \dt 

The connections R = — or u = ef R between the particular solutions of the two equations are 
also basic results of the factorization method together with the construction of the so-called 
factorization partner (or Darboux-transformed) equation of Eq. (T5]) obtained by reverting 
the order of the factorization brackets: 

(j t - RJ + RJ v = = v" - (R 2 - R')v = = v" + (f + 2R')v = . (3) 

Now, we enquire on the implications of a constant shift of the Riccati solution, i.e., 

R s (t) = R(t) + S . 

Rs obeys a Riccati equation of the non-standard form: 

R' s - 2SR S + R 2 s + (f + S 2 ) = J 

which for S = turns into the standard form ([T]). The corresponding linear second-order 
differential equation obtained by substituting Rs = ^| is: 

u" s - 2Su' s + (f+ S 2 )u s = 



and one can immediately see that a particular solution is us = e u. Similarly, v$ = e v is 
a particular solution of the supersymmetric equation 

v'^ - 2Sv' s + (/ + 2R' + S 2 )v s = 

related to the non-standard Riccati equation 

-R' s - 2SR S + R 2 s + {f + 2R' + S 2 ) = . 

However, if we consider the shifted Riccati solution R$ directly in the factorization brackets 
we obtain the following pair of supersymmetric equations: 

&'-[{R + S) 2 + Ii!]$ = , (4) 
[(R + S) 2 -R'}<1> = . (5) 

These equations reduce to the equations for u and v when S — 0. Nevertheless, their 
solutions are not connected to the solutions u and v in the same simple way that us and 
vs are. In the rest of the paper, we will present in full detail the solutions of the latter 
equations in the case of the classical harmonic oscillator for which / is a constant. 

2. The Riccati equation of the classical harmonic oscillator is written as follows: 

R' + R 2 + col = , (6) 

where Uq is the natural frequency of the oscillator. Indeed, using the change of function 
R(t) = — in Eq. (1) we get the harmonic oscillator second order differential equation 

u + UqU = , = const , (7) 

with the general solution u — C\ cos u^t + c 2 sin u^t and the superposition constants c\ and 
C2 fixed through the initial conditions. 

From the two linearly independent particular solutions U\ ~ coscoot and u<i ~ sinu^, 
suppose we choose u\ in the definition of R. Then, one gets 

R Ul = —ujq tan oj t . 

We consider now the factorization partner equation of Eq. (J7|), which is given by 

{li ~ Rui ) {jt + v = ^ v " + ^ v = ' ( 8 ) 

3 



where 

u 2 v {t) = -wjftl + 2 tan 2 w„t) (9) 

denotes the periodic in time 'frequency' associated through the reverted factorization tech- 
nique to the constant frequency of the harmonic oscillator. The linear independent solutions 
v i and v 2 are 



Vi 



V 2 



0J t 1 

— + - sm 2u t 



cos u t u COS u t 

We notice that the u\ and V\ oscillator modes fulfill the relationship 



(10) 



U\V\ = ujq = const , 



but not the same happens with the pair (u 2 ,v 2 ) for which u 2 v 2 = J u\ 2 J u\ is a function 
of time. This is a well-known property of second-order linear differential equations related 
through partner factorizations Q|. 



3. We will now investigate the implications of the constant shifting S of the harmonic Riccati 
solution. The factorization: 



d_ 

dt 



+ Rs 



d 
dt 



R s )U = 



leads to: 



U" + S%(t)U = , tt 2 u (t) 



s 2 s 

— 2 + 2 — tan wot 







UJ 



whereas the partner factorization 



d E 



leads to 



v" + n 2 v (t)v = o , n 2 v (t) 



+ i? 5 v = , 



S 2 S ' 
— r — 2 — tan a> t + 2 tan oj t 



(12) 



(13) 



(14) 



(15) 



Equations ([13]) and (TTo|) define two classes of oscillator equations describing W and V 
oscillators, respectively. Having time dependent frequency parameters £lu,v, these two os- 
cillators are parametric chirps. Moreover, a time periodicity of period T = ^- is embedded 
therein together with periodic tangent singularities. According to Eqs. ( 112]) and ( !T4|) these 
two chirps are connected by the same Riccati solution given by: 



-cuq tan ojot + S . 



(16) 
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Obviously, when S — > we have R$ — > R Ul and the oscillators U and V are very similar to 
the unshifted oscillators u and v in this limit. 



The linear independent solutions U\ and IA 2 have the following form 

U 1 {t)~e- i ^\F 1 (l,-%2----e- 2i "A , (17) 
V ^0 / 

U 2 (t) ~ e^Fi f — , -1; — ; - e -« ) = e ^ st (2 cos 2 u t -ism2u t) ~ e 5t cos^t , (18) 

where = (1 — fj|Vo an d 2-^1 (a, 6; c; z) is the hypergeometric series. Unlike U 2 (t), the 
hypergeometric function in U\(t) cannot be reduced to elementary functions. The simple 
convergence condition dt(a + b — c) < for the hypergeometric series is fulfilled for these 
solutions for all real values of t. We discuss more about the mode Ui(t) in the Appendix. 
Moreover, one can notice immediately that the solutions (jl7p . (jl8j) have a convenient 
Floquet-Bloch form. In addition, the parameter S affects only the period of the phases 
e ±in s t Q f so i u ti ri S but not that of their periodic part. 



On the other hand, the linear independent solutions V are given by: 



e -st 



Vi(t) 7 , (19) 

COS UJ t 



V 2 (t) ~ e st ( U ° + 2Su sin u t + 2S 2 cos u t ) . (20) 
\cos u t J 

Both solutions V\{t) and V 2 {t) are singular at t — and at their odd multiples 

because of the presence of the cosine function in the denominator. Again, one can find a 
pair of solutions (U, V) for which UV = const, in this case (U2, Vi). 

The solutions ([IT]) . ([IB"]) and the exponential factor e st in (fl9l) . (l20l) are bounded if and 
only if the "quasifrequency" Qs has a real value, or equivalently 

wo-iSeR. (21) 

Taking into account that uoq G R, the last condition reads as S = is, s £ R. Notice that 
for a purely imaginary shift parameter, the chirp frequencies and f2y are related through 
&uv(—t) = &u V (t), where * denotes the complex conjugation operation. Such a property 
is called parity-time (PT) symmetry in quantum mechanics, where it has been thoroughly 
investigated since it has been found that the spectra of the PT-symmetric operators are 
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real [lOj. Additionally, by inspecting the solutions (fT7|) , f|T8|) and f|T9|) .f l20|) we note that 
they are periodic for s = (2m — m = 0,±1,±2, ... and antiperiodic for s = 2mojQ, 

m = ±1,±2, ... . 

Plots of the S = 5i periodic and S = Qi antiperiodic solutions are presented in Figures 
1,2 and 3,4, respectively, all of them for u$ = 1. The U. modes are regular in both their real 
and imaginary parts, while the periodic and antiperiodic V modes have periodic singularities 
in their imaginary parts and real parts, respectively. 

For imaginary S, the frequency parameters are: 



n 2 u (t) = ul 

and 



s 2 . s 



1 H 7, + 2i — tan u^t 



CO 



s 2 



2 tan uj t + 2i — tanc^t 



(22) 



(23) 



Equation fT22l) shows that the time-dependent (chirp) feature is present only in the imaginary 
part of the frequency parameter for the U oscillators, whereas their real part looks very 
common. On the other hand, the V oscillators are chirpy in both their real and imaginary 
parts of their frequency parameters and their real part could be even negative. 

3. In summary, in this Letter, we introduced a class of chirp oscillators that are related to 
the classical harmonic oscillators through a constant shift of a Riccati particular solution of 
the harmonic oscillator. Non-trivial results are obtained when the shift parameter is purely 
imaginary and in particular we provide the values of the parameter when the oscillatory 
modes are strictly periodic or antiperiodic. 

We also found that this simple problem displays PT symmetry. Simple classical problems 



with PT symmetry have been also discussed by Bender [10( . However, the way we introduce 
the PT symmetry is by a constant imaginary shift of a particular Riccati solution of the 
problem which is different of the standard procedure in which one employs anharmonic 
imaginary contributions in the potential function. 

One may ask in what particular context these shifted oscillators could be found. Indeed, 
it is not easy to think of a direct physical occurrence of the Riccati solutions. However, we 
have found recently that in barotropic Friedman-Robertson- Walker cosmologies the Hubble 
parameter is just the solution of a special type of Riccati equation, which is the same as 
that of a classical harmonic oscillator unless a linear expression of the adiabatic index of 
the cosmological fluid substitutes the natural frequency ujq of the harmonic oscillator and 
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Fig. 1: The real (solid lines) and imaginary (dotted lines) parts of the periodic solutions U\{t) (in 
red) and 14.2(€) (in blue) for the shift parameter S = 5i. 



instead of the Newtonian time one works in the conformal time Another more terrestrial 
application could be in the context of parametric excitation which results from a periodic 
variation in an energy storage element in a system. To see this, it is sufficient to write the 



shifted oscillator equations for imaginary S in the form 12] 



^ + u 2 U = -ulh{t)U , ^ + s 2 V = -s 2 g(t)V , (24) 

where h(t) = 2ia tanuot — a 2 and g(t) = 2icr~ l tanu^t — 2cr~ 2 tan 2 o> t — < J ~ 2 {p — ^)) respec- 
tively. The equations in (|24p represent the parametric excitation of simple harmonic oscilla- 
tors (or, bandpass filters), of natural frequencies u and s, driven by the signals —u 2 h{t)U 
and — s 2 g(t)V proportional to their responses U and V, respectively. Thus, parametric exci- 
tation experiments with the particular pumpings h and g are a way of studying these types 
of shifted oscillators and their correlations. 
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Fig. 2: The real (solid line) and imaginary (dotted lines) parts of the periodic solutions V\(t) (in 
red) and V2(t) (in blue) for the shift parameter S = hi. 

Appendix: Remarks on the solution Ui{t) = e~ iQst 2 Fi ML, -g; 2 - g; _ e - 2 ^ oi j 



Formula 15.3.4 in Abramowitz and Stegun 13[ (in Linear transformation formulas, page 
559) reads: 

F(a, b- c; z) = (1 - z)~ a F U, c - 6; c; — ^- ) 



Using Eq. (1251) for the HGM series of the U\{b) we get: 

^ ( i. --; 2 - -^ ot ) = ^7 ^ ( i. 2 ; 2 - 7— W ) ■ ( 26 ) 

oj / 2 cos V 1 + e Awo * 



Thus: 

"'W-SisjLs^C 1 ' 2 ' 2 -^!^)- (27) 

On the other hand, we can use the following relationship: 

(1 - z)- p = 2 F 1 (p,&; 6; z) (28) 
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Fig. 3: The real (solid lines) and imaginary (dotted lines) parts of the antiperiodic solutions lA\(t) 
(in red) and U%(f) (in blue) for the shift parameter S = 6i. 

for p — 1 and S = in (11), which leads to: 

U x (t) = - — - — 2 cos u t e~ iuot = e~ iuot (29) 
2 cos uot 

as expected. 



However, if S ^ 0, we get the complicated HGM formula: 

Ul {n) = -^—T ( 2 --)y , (w+1) * - (1 + e 2 ^*)" n • (30) 
U> 2cosa;ot V ^ I> + 2 - g) 1 
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Fig. 4: The real (solid lines) and imaginary (dotted lines) parts of the antiperiodic solutions V\ (t) 
(in red) and V2(t) (in blue) for the shift parameter S = 6i. 
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